Assisted hopping in the Anderson impurity model: A flow equation study 
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We investigate the effect of assisted hopping in the Anderson impurity modei. We use the flow 
equation method, which, by means of unitary transformations, generates a sequence of Hamiltonians 
in order to eliminate the assisted hopping terms. This approach yields a renormalized on-site energy 
e* d , a renormalized correlation energy U* , and other terms, which include pair hopping. For some 
parameter values, the initial Hamiltonian flows towards an attractive Anderson model. We argue 
that this result implies a tendency towards local pairing fluctuations. 
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I. INTRODUCTION. 

The properties of magnetic impurities in metals, and of 
small dots attached to metallic leads share many charac- 
teristics, due to the similar role played by electron elec- 
tron interactions. The electrostatic repulsion leads to 
the formation of local moments^, which are quenched 
at sufficiently low temperatures. Quantum dots are more 
complicated structures^ and an effective Kondo Hamil- 
tonian can only be defined in the limit of large spacing 
between the electronic levels within the dot. Then, the 
main physical processes at low temperatures are due to 
the changes of the occupancy of the level closest to the 
Fermi energy of the leads. Even if we assume this restric- 
tion, terms beyond the Kondo Hamiltonian can arise due 
to the finite extension and inhomogeneities of the elec- 
tron states inside a mesoscopic quantum dot. The sim- 
plest of these additional interactions is an assisted hop- 
ping term, extensively studied in relation to bulk cor- 
related systems 5 * 6 , ft can be shown that, within a sys- 
tematic expansion in the intradot conductance, the lead- 
ing correction to the intradot capacitance, or averaged 
Coulomb interaction, has such a formi The inhomoge- 
neous screening and effective "orthogonality catastrophe 
" which give rise to this term also lead to interesting non 
equilibrium effects in the metallic limit, when the level 
spacing is much less than the temperature^. 

The present work analyzes the model of a correlated 
impurity, or quantum dot, in the limit where the low tem- 
perature physical properties are determined by a single 



electronic state within the dot, using the flow equations 
first proposed in Refs. andllll The method transforms 
the initial Hamiltonian into a family of related Hamiltoni- 
ans by means of a sequence of unitary transformations. If 
adequately chosen, these transformations lead to simpli- 
fied models, in a similar way to a Renormalization Group 
transformation. The method has proven useful for a va- 
riety of problemsi2ii£, and it has also been used to obtain 
the Kondo Hamiltonian from the Anderson modeli^. 

In the following, we will apply the flow equation tech- 
nique to the Anderson Hamiltonian with assisted hop- 
ping. This problem has been analyzed in£ using mean 
field techniques, but there are no more accurate stud- 
ies in the literature. Thus, besides the intrinsic interest 
of the model, this work provides a test case where the 
flow equation method is used to obtain information not 
available by other techniques. 

The next section presents the model, and describes how 
the flow equation method is applied. The results are 
discussed in Section III. Finally, Section IV analyzes the 
main properties of the model, as derived by this method. 



II. THE MODEL AND FLOW EQUATIONS 
A. The model 

We will study the Hamiltonian: 



H — Hk + 'Himp + ^-hyb H - ^assisted 

Hk = J2 e k c k,s c k,s 

k,s 

1-Limp = (?d n d + U°n d ,+n d - 
U hyb = J2Vk(cLds + dtc k , s ) 

k,s 

Uassisted = ^W"n d . s (cl _ s d- s + dL s C k -s) (I) 
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where Ha = T~Lk + 7~timp + T~thyb is the Anderson Hamil- 
tonian, and the assisted hopping terms are included 
in Tiassisted- We have also defined rid,s = d\d s and 
n d = n d ,+ + n d . 



B. The flow equation method 

We now perform a sequence of infinitesimal unitary 
transformations: 



dm = [ v , n] 



(2) 



These transformations are characterized by the anti- 
Hermitian generator 77 = —77^. The parameters of the 
initial Hamiltonians become functions of the flow param- 
eter £ and, if there are no upper indices on the parame- 
ters, it will be implied that they depend on £. 

Our main goal is to reduce the initial Hamiltonian, Eq. 
to another whose physical properties are well under- 
stood. We assume that this is achieved if the assisted 
hopping term is made to flow to zero. Hence, we impose 
the requirement 



lim Ti ass isted(() — > 



(3) 



This condition is satisfied if 77 = [Hk ^assisted]- More 
generally, we will employ the following generator: 



^ Vknd,. 

k . s 



(4) 



This expression differs from the generator used to cancel 
the hybridization with the band, V k in Eq. (JTJ, in that 
it contains an additional operator related to the localized 
orbital, n ds . It is interesting to note that for the analytic 
treatment of the flow equations the explicit expression of 
the parameters rj k is not needed. We only require that 
Eq. © holds. 

The physical properties of the impurity are determined 
by the flow of the on-site energy e d and the correlation 
energy U° as I — > 00. Notice that the one particle hy- 
bridization in the Hamiltonian, Eq. JQ|, is not changed 
during the flow, i.e., 



dm 



hyb 







(5) 



This is a consequence of our choice of the generator given 
in Eq. 

The commutator [77, Ti] generates new interactions 
which have to be included in the Hamiltonian H(£), 
which, in turn, will lead to other interactions. This in- 
finite hierarchy needs to be decoupled in order for the 
method to be useful. This point is carefully discussed 
ini£, and we will follow the same procedure. We do not 
include the new operators in the flow, although they will 
appear in the fixed point Hamiltonian. The role of these 
"marginal" operators will be discussed in section IIIC. 



The resulting flow equations thus take on the simple 
form 



dze d = 2j2vk(V k °+W k )n° k 

k 

d e U = -4j2 7 lk(V k ° + W k ) 



d e W k 



-{e° k -e d -U) m 



(6) 
(7) 
(8) 



where we have introduced the occupation number nS = 
[e^H + which arises from normal ordering the ne- 

glected operators^. The procedure used here explic- 
itly breaks particle-hole symmetry. Note that there 
is no renormalization of the energies of the conduc- 
tion electrons, e°, which is a typical feature of impurity 
problems^. 

From Eq. @, we obtain the general relation between 
the flow of the assisted hopping amplitudes and the gen- 
erator of the infinitesimal transformations, 



Vk 



d e W k 



e d -U 



(9) 



Inserting this expression in Eqs. © and we obtain 



deed 



-2E 



(v k ° + w k )d e w k 



d e U = Aj2 



k c fc 

rO 



e d -U 



(V k ° + W k )diW k 

4 -td-u 



Note that the above equations imply that 



,0 



F(e d 



c d - c d — \ c d T U ) 

U* -U° = G{e" d + U°) 



(10) 
(11) 

(12) 
(13) 

where e* d = e d (£ — 00) and U* = U(£ — 00), and F and G 
denote universal functions. These functions only depend 
on the coupling constants V k ° and W®. 

The above flow equations can also be treated semi- 
analytically by substituting the ^-dependent parameters 
e d and U in the denominator by their fixed point values 
e* d and U* . This strategy has proven to be successful 
in several related modelsiA2ii& We can then integrate the 
equations (fTUfl and ltTT|) from £ — to £ — 00, using the 
constraint that lim^oo W k — > 0. 

We define 



J(e) = (2|F fc °W fc °| - W° k W° k ) S(e - e° fe ) 



(14) 



where we assume that the hybridization V k and the as- 
sisted hopping amplitude W k differ in sign 5 . The spectral 
function is positive for the physically relevant parameter 
regime, where the hybridization amplitude is larger or 
comparable to the assisted hopping amplitude. With the 
above definition we obtain the final result 

J(e) 



de- 



4-u 



:/(e) 



U* - U° = 2 / de 



J(e) 



e-e* d -U* 



(15) 
(16) 
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where we have defined the Fermi function /(e) 

^-114 



III. RESULTS 



(e^ 6 + environment and facilitates the computation. 

To obtain the universal functions F and G from Eqs. 
(|12[1 and (|13|) , we numerically integrate the flow equations 
|0- 0| with 



A. Effective parameters 

In order to make the analysis more quantitative, we 
need to specify the density of states in the surrounding 
medium, described by the function J(e), Eq. (|14|) . The 
width of the conduction band in the leads is 2D, and the 
Fermi energy is ep — 0. We also define 



Vk = (4 ~e d - U)W k 



(18) 



V 



This choice guarantees that Wk vanishes exponentially 
for el ^e* + U* and algebraically for e° = e* + U*M By 
this, we introduce an energy-scale dependent decoupling 
scheme with respect to the renormalized energy of the 
double occupancy of the dot, e* d + U* . 

We now turn to the semi-analytical solution. The spec- 
tral function of Eq. I|14|) reads 



w = JE(^ ) 2 . 



np(e F )(2\V?Wl\ - {W kF f ) = nJ(e F ) 



(17) 



j( ( )Ee(D 2 - f 2 )-^- £ 2 , (19) 



We further assume that both coupling functions, V k ° and with the resonance T = 2(2VW — W 2 )/D. Evaluating 
W®, are semi-elliptic. This choice resembles a realistic Eqs. I|15l) and (|16fl we have to distinguish two cases: 



For \e* d + U*\ < D, we obtain the following self-consistent equations: 



u* - u° 



r 

2D 

2T 
~D 



9 9 l If* -I- 77*1 

(c3 + U*) --D + -Jd 2 - (e* + U*) 2 In ' d 1 = 

Kd ' 7T ttV yd ) d - ^d 2 -{e d + u*y 
(4 + u*) 



For \e* + U*\ > D, we obtain: 

4 + u* 



r 

2D 



-D~sgn(e d + U*)J(e d + U 



D 2 ( 1 arcsin — 

7T (e 



D 



U*) 



u*-u° = 



(e* d + U*) - sgn(eS + U*)^{e* d + U*f 



D 2 



(20) 
(21) 

(22) 
(23) 



In Fig. Q the numerically exact universal functions 
F and G of Eqs. (|12fl and (|13l) are shown as func- 
tions of the scaling variable e d + U° (solid line). The 
coupling strengths were chosen to be V/D = 0.23 and 
W/D — 0.17ii£ The semi- analytical solution (dashed 
line) given in Eqs. (|20l - 123|) qualitatively agrees with 
the exact solution. The main difference is that the kink 
in U* — U° near e d + U° ~ which is present in the 
numerical solution is not adequately reproduced by our 
analytical approach, Eq. JHJ. This is due to the cruder 
decoupling scheme which does not incorportate energy- 
scale separation.lt is interesting to note that at this point, 
the numerical solution yields e* d + U* = 0. It will lead to 
significant differences for the phase diagram which will 
be discussed in the next subsection. 



B. Critical correlation energy 



Eq. lO yields 



U* 



U° - 2e* d T/D 



(24) 



1 + 2T/D 

Hence, for e* d > 0, the renormalized interaction can be- 
come negative. 

We thus define the critical initial correlation energy 
U® such that the renormalized correlation energy U* be- 
comes zero, i.e., 

[/*(t/°| e °)=Lo . (25) 

U® depends on the inital on-site energy e d and we further 
restrict to be positive. 
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FIG. 1: The relative fixed point values t* d — e d (left hand side) 
and U* — (7° (right hand side) as a function of the scaling 
variable e° d + U° for V/D = 0.23, W/D = 0.17 at T = 0. 
The solid lines follow from the numerically exact solution, 
the dashed lines from the semi-analytical approach. 



FIG. 2: Left hand side: The critical correlation energy as 
function of e d following from the numerical (solid line) and 
semi-analytic (dashed line) solution with V/D = 0.23 and 
W/D = 0.17 at T = 0. Right hand side: The critical correla- 
tion energy as function of e d following from the numerical 
solution in the regime where |e^/Z?| < 0.05. 



The critical correlation energy U® is enhanced for pos- 
itive on-site energies e§ > 0. This property can most 
clearly be seen from the analysis of Eqs. (|15f) and (|l(i|> 
for T = oo. 2 ® Then, the Fermi function is constant, i.e. 
/(e) = 1/2, and the integrals yield simple expressions. 
For the semi-elliptical coupling function, we obtain for 
e° d <D + r/2 



2r 



D-r/2 



(26) 



On the left hand side of Fig. |3 the critical corre- 
lation energy U® is shown as a function of the initial 
on-site energy e d as it follows from the numerical (solid 
line) and semi-analytical (dashed line) solution. The cou- 
pling strengths were again chosen to be V/D = 0.23 
and W/D = 0.17. The main difference between the 
two curves stems from the spike of G(e° + U°) around 
e° d + U° rs 0, seen on the right hand side of Fig. 2] 



C. Pairing correlations. 

It is now interesting to consider the possibility of lo- 
cal pairing correlations. Pair formation does not rely on 
particle- hole symmetry, in fact, asymmetries can actually 
enhance this phenomenon^ But local pairing can only 
occur if - in addition to an effective negative correlation 
energy U* < - the renormalized on-site energy e* d lies 
sufficiently close to the chemical potential /i = The 
local pairing regime is thus determined by the tradeoff 
between a high critical correlation energy £7? and an ef- 
fective on-site energy e* d sufficiently close to the chemical 
potential. 



To discuss the phase diagram of this local pairing 
regime we will only consider the numerical solution. This 
approach is more reliable since it incorporates a system- 
atic decoupling of the assisted hopping term, starting 
with large and ending with small energy differences with 
respect to e* d + U* . The results are presented on the 
right hand side of Fig. |2 where the critical correlation 
energy U® is shown as function of e d in the regime where 
< 0.05. Furthermore, we distinguish between pos- 
itive (+) and negative (A) renormalized on-site energy 
f* 

Our analysis predicts local pairing for values of e d 
which favor double occupancy, and £/?<0.1. The mean 
field analysis in Ref. gives local pairing for t/?<0.05 
and a similar range of values of e d (note that we are us- 
ing here comparable hybridization with the band). 

It is interesting to note that the strong renormaliza- 
tion of U* for e d + U° s» only occurs at large values 
of £ 3> l/D 2 . Therefore, it is not present in the semi- 
analytical solution which does not incorporate energy- 
scale separation. Furthermore, there are fixed points 
for which e d ~ U* « 0, i.e., a degeneracy between the 
zero-charge and two-charge state of the impurity. Since 
the marginal interactions also include pseudo spin-pseudo 
spin couplings, this might give rise to the Kondo effect 
in the charge channel, i.e., in the mixed valence regime. 



D. Marginal interactions. 

We now analyze the marginal terms, which have been 
neglected so far. One term will be interpreted as a renor- 
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malization of the hybridization term V k °. This identi- 
fication is crucial since the assisted hopping term W k 
contains an effective hybridization term which must be 
present in the fixed point Hamiltonian. Another term 
will explicitly favor local pairing. 

As in the Anderson model, the flow generates spin-spin 
and pseudo spin-pseudo spin couplings: 



T~tsw - - E V k*k> 



k,k> 



^t^ipk' ■ +h.c. 



(27) 



In addition to this interaction, the following term is 
also generated: 

K lr r = \Y. V W ( : ^l*>Pk> ■ +h.c?j + ~e* d {n d - 1) (28) 
k,k' 

The existence of these terms is independent of the as- 
sisted hopping interaction and has been discussed in de- 
tail in relation with the Anderson models. There, Hsw 
gives rise to the well-known Kondo behavior in the spin 
and charge sector, respectively. 

Ti.i rr is usually neglected in the fixed point Hamiltonian 
since the first term of Eq. (|28l) represents potential scat- 
tering of the conduction electrons and will be irrelevant 
for our discussion. The last term renormalizes the on-site 
energy by e* d = J^k^Ck/^- For the semi-elliptical cou- 
pling function l(T§)> . we obtain e* d /e* d » T/2D according 
to Eq. (|3U|) within the semi-analytical treatment. Our 
approximation - not to include newly generated coupling 
terms into the flow of the Hamiltonian - is thus consistent 
for T/D <C 1. As in Ref. 1 1 2t we will ignore all terms of 
Eq. I|28|) in the fixed point Hamiltonian and assume that 
they can be taking into account by redefining the initial 
parameters. 

Finally, the commutator [77, Tiassisted] also gives rise to 
the following pairing term: 



Impair — 



E 



wt k , 



c k,s c k' 



ed- s d s + h.c. 



(29) 



which has the same shape as some of the contributions 
in Hsw, Eq. flT7)> . 

The additional contributions of Eq. (|27l - 129|) need to 
be included in the flow equation scheme. As these terms 
contain an even number of operators associated to the 
localized level, but the generator, Eq. |0J, contains an 
odd number, these terms cannot directly change the flow 
equations of the parameters U and Cd, Eas. l|6l7|) . They 
can, however, change the flow of the assisted hopping 
term, Eq. (jSJ), which will eventually modify Eqs.© and 
J7J). We will neglect this second order modification of 
the flow equations for U and e^. This approximation was 
initially introduced ini£ and is based on the fact that 
V kk , = and W k ,k' = for t = 0. Note that, provided 
that the flow equations are unchanged, the presence of 



these "marginal" terms in the fixed point Hamiltonian 
does not alter the main conclusions of the paper. 

We now approximate the magnitude of these terms 
at the fixed point, neglecting their influence on the flow 
equations. We obtain: 



d£r) k {V k °, + W k ,] 



2|W«,| 



u* 



d£r] k W k 



w k w k' 



2e k < 



U* 



(30) 
(31) 



As mentioned earlier, Hsw arises in the analysis of 
the Anderson Hamiltonian, due to the effects of the hy- 
bridization between the state at the impurity and the 
conduction band, and can be derived by means of a 
Schrieffer- Wolff transformations, . Hence, we can inter- 
pret these fixed point interactions as arising from an ef- 
fective hopping term in the initial Hamiltonian, which, 
indeed, can be obtained from a mean field decoupling of 
the assisted hopping term. TL pa ir will enhance the ten- 
dency towards pair formation. 



IV. CONCLUSIONS. 

We have studied a model of an impurity, or a quan- 
tum dot, with a single electronic level, an on-site repul- 
sion term, hybridization, and assisted hopping. The flow 
equation method allows us to map the model onto an 
effective Anderson model without assisted hopping, and 
with additional Kondo couplings and pair hopping terms. 

The most interesting outcome of our calculation is that 
the value of the on-site electron-electron interaction can 
become attractive, because of the effect on the flow of 
the initial assisted hopping term. Although the Ander- 
son model with attractive interactions has not been as 
extensively studiecU&i^ as the repulsive Anderson model, 
we think that our mapping suffices to give a qualitative 
description of the physical properties of the model for 
the interesting case of a relatively strong assisted hop- 
ping term. Note that, in addition to a renormalized on- 
site interaction term, terms which describe pair hopping 
between the impurity and the conduction band are gener- 
ated. We expect that significant local pairing correlations 
will develop in this regime. 

This result is in agreement with general arguments 
which suggest that assisted hopping interactions tend 
to favor superconducting ground states*, and with mean 
field studies of the same impurity models. A related ten- 
dency towards phase rigidity can be found in the limit 
where the level spacing within the impurity is negligible 
with respect to the other parameters^. 

The model does not have electron-hole symmetry, and 
local pairing is mainly favored when the number of elec- 
trons at the impurity exceeds one (half filling) , and when 
the direct hopping and the assisted hoping terms have 
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